Reciprocity and the scattering matrix of waveguide modes 



Guro K. Svendsen, 1 Magnus W. Haakestad, 2 and Johannes Skaar 1 

1 Department of Electronics and Telecommunications, 
Norwegian University of Science and Technology, NO-7491 Trondheim, Norway 
2 Norwegian Defence Research Establishment (FFI), NO-2027 Kjeller, Norway 
(Dated: January 14, 2013) 

The implications of the Lorentz reciprocity theorem for a scatterer connected to waveguides with 
arbitrary modes, including degenerate, evanescent, and complex modes, are discussed. In general 
it turns out that a matrix CS is symmetric, where C is the matrix of generalized orthogonality 
coefficients, and S is the scattering matrix. Examples are given, including a scatterer surrounded 
by waveguides or free space, and discontinuities of waveguides. 



I. INTRODUCTION 

The Lorentz reciprocity theorem is a famous and very 
useful result in electromagnetics [THS]. In particular, it 
has important consequences in scattering theory, relat- 
ing the transmission coefficients when the position of 
the source and detector are interchanged. For orthog- 
onal, normalized, propagating modes in waveguides or 
free space, the forward and backward transmission coef- 
ficients due to a scatterer with symmetric permittivity 
and permeability tensors, are equal. This symmetry has 
proved to be useful in a number of situations in physics 
and electrical engineering. For example, the symmetry 
implies that the transmission coefficients through win- 
dows, thin-film filters, or fiber Bragg gratings are equal 
from opposite sides. The symmetry property also shows 
that one-way components such as isolators and circula- 
tors, cannot be made of reciprocal materials. 

With the demand for device miniaturization, coupling 
to and from evanescent modes becomes increasingly im- 
portant. For instance, when considering the end-facet 
reflection and transmission of nanowire lasers, or the re- 
flection and transmission of nanowire Bragg reflectors, 
the evanescent waveguide modes play an important role 
[6l [7] . For certain metal-dielectric waveguides, there are 
modes with complex propagation constants 8, 9], even 
though the material is lossless. In the presence of degen- 
eracy, modes may not be orthogonal. With the recent 
interest in plasmonic, metamaterial, and photonic crys- 
tal waveguides it is thus of interest to see if there are 
symmetry relations for the scattering matrix when com- 
plex and/or degenerate waveguide modes are included. 

While symmetry relations are well known in the case 
of propagating modes, the case with evanescent modes is 
less developed in the literature. Vassallo discussed the 
symmetry relation for the scattering matrix in the pres- 
ence of evanescent waveguide modes g]. For the special 
case with free space on each side of the scatterer, Carmi- 
nati et al. found the reciprocity relation for the scatter- 
ing matrix associated with plane, electromagnetic waves 
5J. Earlier treatments did not include the possibility of 
complex modes, and nonorthogonal modes. 

Here we will provide a self-contained treatment of the 
Lorentz reciprocity theorem and its implication for the 



scattering matrix or transmission coefficients between 
propagating, evanescent, and complex modes. We also 
consider the effect of degeneracy in detail. The scatterer 
is contained in a black box with one or more waveguide 
inputs, and/or free space on one or both sides. As stated 
in detail in the next section, the analysis is quite general; 
it only requires the usual assumptions for the Lorentz' 
reciprocity theorem to be valid for the scatterer, and the 
additional assumption that the input / output waveguides 
are described by scalar permittivities and permeabilities. 
Examples are given, which demonstrate the reciprocity 
relations for a discontinuity between two waveguides. 

It will turn out that it is not the scattering matrix S, 
but the product CS, that is symmetric. Here C is the ma- 
trix of generalized orthogonality coefficients Q . This re- 
sult is useful when designing or analyzing structures with 
nontrivial waveguides, including plasmonic or nanowire 
waveguides. For example, when designing waveguides 
with Bragg gratings or other inhomogeneities along the 
waveguide axis, it is useful to relate the reflection and 
transmission coefficients from the two sides, in order to 
identify transfer matrices of the different sections. Also, 
it is important to have the symmetry in mind, to be aware 
of the inherent limitations for structures satisfying the 
conditions of the Lorentz' reciprocity theorem. 



II. LORENTZ RECIPROCITY AND THE 
SCATTERING MATRIX 

Consider a volume v with a linear, time- invariant, and 
spatially nondispersive medium. The medium may be 
lossy, temporally dispersive, anisotropic, and inhomoge- 
neous. Let e and fi be the permittivity and permeability 
tensors in v, respectively. We assume that e and /i are 
symmetric; this property follows from thermodynamic 
relations under very general conditions [10 . The con- 
ductivity, if any, is included in the (possibly complex) 
permittivity e. The volume v is seen as a black box with 
an unknown system; however in addition to the assump- 
tions above it is known that there are no current sources 
in v. Let {E 1 , H 1 } and {E 2 , H 2 } be the electromagnetic 
fields in v, resulting from excitations 1 and 2, respec- 
tively, whose sources are located outside v. The Lorentz 
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FIG. 1. The "black box", or scatterer, is contained in a vol- 
ume v with boundary s. (a) The general case with several 
waveguides connected to the surface s. (b) The special case 
with two waveguides, one on each side of v. The surface s is 
regarded as two infinite planes. 

reciprocity theorem reads (see Appendix A): 

jf(E* x H 2 ) -dS = j(E? x H 1 ) -dS, (1) 

where the closed surface s is the boundary of v. 

We now consider a situation where the volume v has 
several input /output waveguides (Fig. 1). Each of the 
waveguides is described by a permittivity and permeabil- 
ity which are uniform along the waveguide's axis. The 
waveguides may be lossless or lossy. We assume that any 
solution to Maxwell's equations in the waveguides are 
expressible as a superposition of modes: Bound modes, 
radiation modes, and/or evanescent and complex modes. 
The waveguides are assumed to be isolated, in the sense 
that their modal fields do not overlap. 

An excitation in one or more waveguides may lead to 
forward and backward propagating waves in each mode. 
Here, "forward" refers to the input direction, while "back- 
ward" refers to the output direction. A mode is referred 
to as "forward propagating" if its power flow is directed 
in the forward direction (for propagating modes) or if it 
decays in the "forward" direction (evanescent or complex 
modes) , see Appendix [B] Let the modal amplitudes of 
the forward and backward propagating waves be dj and 
bj , respectively. This set includes all modes in all waveg- 
uides. The modes are labeled by j, which is discrete in the 
discrete part of the mode spectra, and continuous other- 



wise. Alternatively, we can ensure that the entire mode 
spectra are discrete by imposing artificial boundary con- 
ditions, such as periodic or metallic boundary conditions. 
For example, if a dielectric waveguide is surrounded by a 
metallic box, the continuous set of radiation modes gets 
discretized. 

When two or more waveguides do not have parallel axes 
(Fig. 1(a)), there is a fundamental problem with the ra- 
diation modes associated with dielectric waveguides. In 
principle, such modes extend infinitely in the transver- 
sal waveguide direction, and will therefore overlap spa- 
tially, as opposed to our earlier assumption about isolated 
waveguides. There are at least two remedies. The first 
remedy is to consider supermodes rather than modes of 
separate waveguides. In other words, one considers the 
true modes of the composite structure consisting of two 
or more waveguides. The second, perhaps more practical 
remedy is again to use artificial boundary conditions, lim- 
iting the transversal dimensions of the modal fields. In 
the following we will use notation for the case where the 
modal fields have finite transversal dimensions; the case 
where the modal fields extend infinitely in the transversal 
direction is described by substituting sums by integrals, 
and Kronecker deltas by Dirac delta functions. 

At the surface s we can write the transversal compo- 
nents of the electric and magnetic fields as follows: 

E'=EM + ^), (2a) 

j 

h 1 = EM5-^). (2b) 

Here {e^ , hj} denotes the modal fields of mode j, and the 
superscript t stands for the transversal component. The 
sums run over all modes in all waveguides. Note the mi- 
nus sign of the backward propagating magnetic field; this 
is due to the usual mode convention in the literature (see 
Appendix [b]). In ^ the spatial dependence exp(±if3jZj), 
where /3j is the propagation constant of mode j, and Zj 
is the associated waveguide axis, has been absorbed into 
the modal amplitudes ctj and bj. 

Since the system in v is linear, we can express the out- 
put propagating modes using a scattering matrix Sjk- 

bj=^2S jk a k . (3) 

k 

Defining 

J ej xh fe -dS = C ife , (4) 

the properties of the matrix Cjk are found to be as fol- 
lows: When j and k refer to modes in different waveg- 
uides, orthogonality Cjk — is trivially fulfilled since 
the associated modal fields are assumed to be non- 
overlapping. For the nontrivial case with j and k in 
the same waveguide, it is proved in Appendix [B] that 
Cjk is symmetric, and that Cjk = for modes j and 
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k with distinct propagation constants. Thus Cjk is a 
block-diagonal, symmetric matrix. We here note a gen- 
eralization compared to earlier treatments in the liter- 
ature: Due to possible nonorthogonality of degenerate 
modes, we permit Cjk to be nondiagonal. 

We now consider two different excitations and their 
associated fields on s. The first has a j = 1 for a certain 
j, putting all other aj's to zero; and the other has a\ = 1 
for a certain k, putting all other a^'s to zero. Here the 
superscripts 1 and 2 denote the two different excitations. 
This means that 



We note that (11) is valid even though there are 



E 



H 



i 



and 



E 2t =et 



i 



(5a) 
(5b) 

(6a) 
(6b) 



Substituting ^ and ^ in the Lorentz reciprocity theo- 
rem 0, and using Q yields 

_ ClmSljSmk + ^ (ClkSlj — CjlSlk) (7) 
l,rn I 

= Ckj — 2_j CmlSljSmk + (CljSlk — CklSlj) . 

Using the symmetry of Cjk we conclude that 

^ CjiSik — ^ CkiSij, (8) 

for all j and fc, i.e., the matrix CS is symmetric: 

CS={CS) T . (9) 

In some cases of interest, we may have orthogonalized 
part of the mode set, while we do not care about the re- 
maining part. This may for instance be the case when 
considering certain bound modes in dielectric waveg- 
uides, or when some of the input /output waveguides 
are irrelevant to the problem under investigation. For j 
and/or k in this restricted, orthogonal set, we can write 
(Appendix [B]): 



x h k -dS = CjSjk, 



(10) 



where the symbol Sjk stands for the Kronecker delta, 
and Cj is a normalization constant. In other words, 
Cjk = CjSjk for j and/or k in the orthogonal set. By 
substitution into ^ we get the important result: When 
mode j and mode k are orthogonal to each other and to 
all other modes, 



nonorthogonal modes in the remaining part of the modal 
spectrum. 

One common normalization choice is to let 



1, for propagating modes, 
±i for evanescent modes, 



(12) 



while modes with complex or zero propagation constant 
may have general complex Cj. Then, when modes j and 
k are orthogonal to each other and to all other modes, 

Sjk = Skj (13) 

for propagating modes, while 

Sjk = ±iS k j (14) 

if j is propagating and mode k is evanescent. 

We now consider some common special cases of the 
general reciprocity relation 

III. EXAMPLES 

A. Orthogonal waveguide modes on the left and 
free space on the right 

For the special case where one of the involved waveg- 
uides vanishes, i.e., we have free space on the right-hand 
side of the volume v in Fig. [ijb) , we can take the modal 
fields on the right-hand side to be complex exponentials, 
or plane waves: 



ej exjp(ik x x + ik y y), 
hj exp(ik x x + ik y y). 



(15a) 
(15b) 



Here the constant vectors and hj express the polariza- 
tion. Using Maxwell's equations we find the connection 
between these constant vectors, and also the dispersion 
relation 



P 2 



e/j,- 



k y- 



(16) 



Here e and fi denote the relative permittivity and per- 
meability, respectively, and w is the angular frequency. 
The modes are labeled by k x , k y , and the polarization 
pol = TE or TM, so when j and k refer to modes in 
free space, they can be seen as a collection of indices, 
j — s- (k x , k y , pol) and k 



we let — j — > ( 



-k x , — k y , pol) and 



(K,k' y , P oY) 

k - 



For convenience 
(-^-fc^pol')- 

We take h^- to be real. Although not assumed here, for 

real e and fi this means that is real for propagating 
modes, and imaginary for evanescent modes. Using an 
appropriate mode normalization, we find the orthogonal- 
ity relation 



e(fc,,fc yiP oi) x h ( _ feii _ fe / pol /) ■ 



dS 



CjSjk — CkSkj- 



(11) 



c(k z ,pol)S(k x - k x )S(ky - k' y )S po y >pol , (17) 
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by direct calculation. Here the integration extends over 
the entire xy-plane, and 



2 IM 



c(fc„pol) = ^Mi 
0, 



for pol = TE, 
for pol = TM, 
for k z = 0. 



(18) 



Taking the other waveguides into account, and assum- 
ing orthogonal modes in them, we obtain Q, with 



We find that (21c) remains valid for the reflection coef- 
ficients on each side when the modal fields are taken to 
be complex exponentials. If there are different homoge- 
neous media on each side, the appropriate electromag- 
netic parameters must be used in (16) and (22). For the 
transmission coefficients we obtain 

c(fe^,pol)S(ieft,fc a .,fc !( ,p l), (right, k' x ,k' , pol') 

= c(fc^,pol )5( ri ght,-fe;,-fe' polO.CLeft.-fca.-fcy.pol)) (24) 



C 



jk 



c(k z , pol)<Jj(_fe), j and k refer to free space modes,where "left" and "right" specify whether the particular 
CjSj k , j and k refer to waveguide modespode is located on the left-hand side or right-hand side 

otherwise °^ ^ e sca ^ ereT - This special case was treated previously 

in Refs. [5JHI]. 



(19) 

The short hand notation <5j(_fc) for j = (k x , k y , pol) and 
k = (k' x , k' , pol') means 



= 5(k x + k' x )S(k y + k' y )5 po i tP0 i' 



(20) 



Due to the nonorthogonality of the degenerate modes 
{k x , k y , pol) and (— k x , —k y , pol), C is not diagonal in this 
case, so we invoke the general version Q of the symme- 
try relation. This gives (here we note that the sum is 
replaced by an integral for the continuous spectrum of 
free space modes): 



(21a) 



c j%(fc*,fcjMPol) = c (k z ,Vo\)S(_ kx _ ky:PO i)j, (21b) 

c(fc z ,P0l)5'(fc :c ,fc !/ ,pol),(fc! c ,fc;, P ol') 

= c(k' z ,pol')S(_ K -fe^poi'W-fe* ,-jfe vl poi), (21c) 



In (21) j and k denote waveguide modes, while 
(k x , k y , pol) and (k' x , k' y , pol') denote free-space modes on 
the right-hand side of the scatterer. Moreover 



efiu> 



>,'2 



','2 



(22) 



Instead of using the complex exponentials ( 15 ) as 



the modes, we can form linear combinations of modes 
with opposite signs of k x , and opposite signs of k y . 
This gives modal fields with components of the form 
cos(k x x) cos(k y y) , cos(k x x) sm(k y y), sm(k x x) cos(k y y), 
or sin(k x x) ain(k y y). It is straightforward to verify that 
the resulting modes are orthogonal, i.e., C now becomes 
diagonal. Thus in this case we obtain 

CjS'jfc = Cf-Si-j, (23) 

rather than (21 ). Here j and k may denote waveguide or 
free space modes. For free space modes cj — c(fc z ,pol) 
and Cfc = c(fe^,pol). 



B. Free space on both sides 

When there is free space on both sides of the scatterer, 
the analysis is similar to that in the previous subsection. 



If the sin/cos type modes are used instead, (23) re- 
mains valid. 



C. Discontinuity between two waveguides 

Consider the discontinuity between two waveguides, 
see, for example, Fig. [2| The goal is to calculate the 
reflection coefficients and transmission coefficients from 
either side, demonstrating that they obey the reciprocity 
relation (11). For completeness and notational consis- 



tency, we include a derivation based on Refs. [H [7] . 

First we let the two waveguides be arbitrary. How- 
ever, for simplicity we use sum notation for the superpo- 
sition of modes, indicating that any continuous part of 
the mode spectrum is discretized using artificial bound- 
ary conditions. Furthermore we assume that the mode 
sets are orthogonal, i.e. they satisfy 



i) x hi 



dS 



[ e 2 1 xh 2 k -dS = c 2 ] 5 lk , 

J A 



(25a) 
(25b) 



where A is the boundary between the two waveguides, 
and cj' 2 denotes the normalization constant of mode j in 



waveguide 1,2 (see (B16)). Define the overlaps 



/ e]xh 2 -dS, 

J A 



* 2 fe = ^e 2 fe xhj-dS. 



(26a) 
(26b) 



We consider the situation where mode j is incident 
from waveguide 1 to waveguide 2, and impose continu- 
ity of the transversal electric and magnetic fields. This 
amounts to solving the equations 



(27a) 
(27b) 
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where rL is the reflection coefficient from mode j to mode 
k, as seen from waveguide 1, and tjy is the transmis- 
sion coefficient from mode j in waveguide f to mode k in 
waveguide 2. Using (25) and poj) we can rewrite p7|: 



c 1 r 1 ■ = ViL* 1 u 

mj / j kj m k ' 



(28a) 
(28b) 



(1) 



Summing these equations we find 



(29) 



This matrix equation can be inverted, yielding a matrix 
Fresnel equation expressing the transmission coefficients 
t\- from the modal fields. 

Repeating the derivation above, we find the transmis- 
sion coefficients from waveguide 2 to waveguide 1 by solv- 
ing 



(30) 



Comparing (29) and (30) as matrix equations, we realize 
that 



ly.2 _ 2.1 



(31) 



in agreement with the reciprocity relation (If). 



The reflection matrix can be found by eliminating the 



transmission coefficients from ( 28 ) , or alternatively, by 



formulating (27) as follows 



$i + X^rl Vl - c 2 t 1 ■ 

J m / * kj km ^m mj s 
k 


(32a) 


* 2 - Vr! - c 2 t 1 

J m / j kj km m mj ' 

k 


(32b) 


In matrix notation, we define R 1 = [rLj and ^> 1,2 = 

[^J. Eqs. (32) can then be combined to the matrix 
Fresnel equation 


R 1T = (* 2 - ^r 1 )^ 1 + ^r 2 )" 1 


(33) 


Similarly we find the reflection matrix 
waveguide 2: 


as seen from 


R 2 = (* 1 + * 2 )" 1 (* 1 -* 2 ). 


(34) 



To demonstrate that these expressions are consistent 



with reciprocity (If), we specialize to two planar, non- 
magnetic waveguides, and an incident TE mode. Both 
waveguides are oriented such that the structure is uni- 
form in the y-direction. Since there will be no coupling 
between TE and TM modes, we can disregard TM modes 
from the analysis. By straightforward calculation we find 



(35a) 
(35b) 



FIG. 2. The discontinuity between parallel plate, metallic 
waveguides. On the left-hand side (I) there is a waveguide of 
thickness d 1 inserted into another waveguide of thickness d 2 . 
On the right-hand side (2), only the waveguide of thickness 
d 2 is present. 



where (3j (/3|) is the propagation constant of mode j (k) 
in waveguide 1 (2), and 



f 



efdx. 



(36) 



cross section 



From ( 36 ) it is apparent that even and odd modes do not 
couple to each other, so we limit the discussion to even 
modes. 

First we consider the situation where the waveguides 
are parallel plate, metallic waveguides, with thicknesses 
d 1 and d 2 , respectively, and common axis z, see Fig. [2] 
The waveguides are filled with vacuum. The modes are 
easily found analytically, yielding analytic expressions for 
(3j, /3 2 , and fyjk- Using (33) and (35), the resulting ma- 



trix c\r\j is calculated and shown in Fig. p| We observe 
that the matrix is symmetric, as predicted T5y reciprocity: 
Eq. ( II ) can in this case be written 



c j r jk 



-k'ky 



(37) 



1 for all k, which 



Using the normalization (12), |cfc| 
means that \rkj\ becomes symmetric. 

Second we consider the end facet of a thin, dielectric, 
planar waveguide of thickness d. Using p3| , (35), and 
( 36 1 we compute the reflection matrix rL . The result- 
ing matrix |cfcrLj is depicted in Fig. |4j and is clearly 
symmetric. In Fig. p a) , we also show the phase of r jh , 
which demonstrates the asymmetry of the scattering ma- 
trix due to presence of evanescent modes. For compari- 
son, Fig. [5](b) shows the phase of the symmetric matrix 

c k r jy- 



IV. CONCLUSION 

We have considered the implications of the Lorentz 
reciprocity theorem for the scattering matrix. The modes 
to and from the scatterer are described as general waveg- 
uide modes, including propagating, evanescent, and com- 
plex modes. The modes may also be degenerate and 
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FIG. 3. The matrix |cfc»"jy[. The numbers used are d 1 = d, 
d 2 = 2d, and tod/c — 20. In the plot, mode 1, 2, . . . , 15 cor- 
respond to even TE modes in the central waveguide on the 
left-hand side, while 16, . . . , 30 correspond to even TE modes 
in the upper and lower waveguides. Modes 1, 2, 3, 16, 17, and 
18 are propagating, while the remaining modes are evanes- 
cent. The plot was generated using 100 modes on each side in 
the calculations, which leads to an asymmetry of 10 -3 for the 
matrix. This error is reduced to 10 -5 by using 1000 modes. 

nonorthogonal. Despite this generality, a simple reci- 
procity relation ([9| can be derived, stating that the ma- 
trix product of the orthogonality matrix C (with elements 
defined by Q) and the scattering matrix S, is symmetric. 
This symmetry relation is expected to be useful in pho- 
tonics and microwave engineering, when designing struc- 
tures involving complex waveguides. 

Appendix A: Lorentz' reciprocity theorem 

Consider a volume v with a linear, time-invariant, and 
spatially nondispersive medium. The conductivity, if any, 
is included in the (possibly complex) permittivity tensor 
e. We assume that there are no current sources in v, and 
that the tensors e and \x are symmetric. Let {E 1 ,!! 1 } and 
{E 2 ,H 2 } be the electromagnetic fields in v, resulting from 
excitations 1 and 2, respectively. From Maxwell's equa- 
tions we have (assuming time dependence exp(— iujt)): 

V • (E 1 x H 2 ) = -E 1 • V x H 2 + H 2 • V x E 1 

= icjE 1 • (eE 2 ) + iwH 2 • (pU 1 ) (Al) 

and 

V • (E 2 x H 1 ) = iwE 2 • (eE 1 ) + icuH 1 ■ (^H 2 ), (A2) 
Using the symmetry of e, we have 

E 1 • (eE 2 ) = J2 Eie kl Ef = J2 E?e lk El = E 2 . (eE 1 ). 

fc,J=l k,l=l 

(A3) 




FIG. 4. The matrix le^r^l. The refractive index of the dielec- 
tric waveguide core is set to 1.5, and uid/c — 1. The modes 
are computed numerically using periodic boundary conditions 
with period L = lOOd. To reduce the effect of the artificial 
boundary conditions, we introduce a small loss everywhere, 
such that e — > e + i0.05. The coefficients for the 50 lowest 
order TE waveguide modes are displayed in the figure. Only 
the first waveguide mode is bound, modes 2, . . . , 31 are un- 
bound and propagating (i.e., with small imaginary parts of 
the propagation constants), and the remaining modes are un- 
bound and evanescent (i.e., with large imaginary parts of the 
propagation constants). The plot was generated using 100 
modes on the waveguide side. On the free space side, for 
simplicity we used complex exponentials as the modal fields. 
Since these fields are neither even nor odd, we included twice 
as many modes there, corresponding to positive and negative 
transverse wavenumbers. When the matrix Fresnel equation 
( |33[ ) is overdetermined, the matrix inverse should be replaced 
by the Moore-Penrose pseudoinverse. Using 100 waveguide 
modes leads to a matrix asymmetry of 10 -3 . Again, the error 
is reduced to 10 -5 when using 1000 waveguide modes. 

Similarly, 

H^O/H 2 ) -H 2 -^ 1 ). (A4) 

We therefore obtain the usual Lorentz reciprocity theo- 
rem: 

V^E 1 x H 2 ) = V-(E 2 x H 1 ), (A5) 
or, in integral form: 

^E 1 x H 2 ) -dS = j((E 2 x H 1 ) -dS, (A6) 

where the closed surface s is the boundary of v. 

Appendix B: Waveguide modes 

For structures that are homogeneous in the z-direction, 
the solutions to Maxwell's equations are expressible as 
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tains 






3 


e* = 


- — z x (Bh* + i^h* 




2 




— z x f^ + iVV) , 
a;// 




1 


e z = 


— z-V* x h\ 
we 







h z = 


-^z-V 4 x e*. 



40 50 




-0.5 



-1.5 



10 20 30 40 50 



FIG. 5. (a) Phase of rL, and (b) phase of c^rL, for the 
dielectric waveguide in Fig. [4] The phase plots clearly show 
that multiplication by the Ck coefficients is necessary to obtain 
a symmetric matrix. 



superpositions of waveguide modes, that is, solutions in 
the form 



E = a(e t + e z z) exp(i/3z), 
H = a(h' + h z z)exp(i/3z). 



(Bla) 
(Bib) 



Here the superscripts t and z stand for the vector com- 
ponents perpendicular and parallel to the z-axis, respec- 
tively, and a is a constant. Assuming a (possibly lossy) 
medium with permittivity e and permeability /i, and sub- 
stituting into the source-free Maxwell equations, one ob- 



UJfl 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



It is therefore consistent to choose "backward propagat- 
ing modes" as related to "forward propagating modes" as 
follows: Under the transformation 8 — > —3, we require 
e* -> e\ e z -> -e z , h z -> h z , and h* -> -h*. 



By inser ting Eq. ( |B2d[ ) in Eq. ( |B2a[ ), and Eq. < \B2c\ in 

Eq. (B2b|, and rearranging the terms, Maxwell's equa- 
tions take the form of a generalized eigenproblem [121 113] 

A\9)=p6\V), (B3) 

where 



A = 



-V* x iV*x 






cj/i - iv* x iy* 



B = 



-zx 
zx 



and 



(B4) 



(B5) 



(B6) 



The inner product between two vectors and \^k) is 
defined as 



[(ej-r-e^+^r-hyds, (b?) 



where the integration area A either is the entire xy-plane 
or the computational area associated with, for example, 
periodic or metallic boundary conditions. It is straight- 
forward to prove that the operator B is hermitian and 
real, and therefore symmetric. Also B^ 1 = B. For gen- 
eral lossy media, A is not hermitian; however we can 
prove that it is symmetric (i.e., A*^ — A) by verify- 
ing the equation (^j\A\^ k ) = ( J 4*|^))t|* fc ). Indeed, 
V • (u x v) = (V x u) • v — u • (V x v), and therefore 



/ (V x u) • vdS* = f (V x v) • udS 

J A J A 



(B8) 



whenever the line integral § gA {u x v) • ndZ = 0. Here fi 
is the outward-pointing unit normal vector to the bound- 
ary dA of A. This condition is satisfied for metallic or 
periodic boundary conditions, and can be justified in gen- 
eral by enclosing the waveguide by a conducting cylinder 
whose radius approaches infinity. 



Transposing (B3), and using the symmetry of A and 
-B, we obtain 



(%*\A = /3{y*\B. 



(B9) 



By considering (B3) for a waveguide mode |^fc), taking 
the inner product with I 1 ?*), we obtain 



(Pj-0 k )(9* j \B\* k ) = O. 



Thus if f3j 7^ /3fc we have 



k 



With the help of (B7), Eq. (Bill can be written 



1 / (ej x h fe + e fe x hj) ■ zdS = 

2 J A 



(BIO) 



(Bll) 



(B12) 



for j3j 7^ Using the transformation properties of 
and hfc as — > ~Pk, we obtain 

- / (-ej x h fc + e fe x hj) ■ zdS = (B13) 
2 J A 

for ;Sj 7^ —pk- Thus, for j3j 7^ ±/3fc we obtain 



e, x h k ■ zdS = 0, 



(B14) 



by subtracting (B13| from (B12). Note that the or- 



thogonality (B14) is valid for all nondegenerate waveg- 



uide modes, including complex modes. For degenerate 
modes, however, we are not guaranteed that (B14| is 



fulfilled, even though the mode set may in some cases 



be orthogonalized, for example using the Gram-Schmidt 
procedure^] In general we therefore write 



x hi, ■ zdS = C 



(B15) 



where Cjk is a block diagonal matrix such that Cjk = 



for nondegenerate modes (J3j 7^ ±/3fc). From (B13l we 
find that Cjk = Ckj provided j3j 7^ —/3k- The latter con- 
dition can always be assumed fulfilled by restricting the 
mode set to forward propagating modes, i.e., by mak- 
ing a mode set consisting of the half of modes satisfying 
(3j 7^ — For gainless waveguides this amounts to pick- 
ing the modes with power flow in +z-direction (prop- 
agating modes), or Im/3 > (evanescent and complex 
modes). For modes with real f3 and zero power flow in 
the z-direction, the forward propagating subset can be 
found by considering a lossy waveguide in the lossless 
limit, or by invoking causality. 

Provided the degenerate modes can be orthogonalized, 
we can write 



ej x hfc 



idS = Cj5 



■ 3 o 3k - 



(B16) 



In principle we can take Cj = 1 for all j by a suitable 
normalization. It is however conventional to let 



1, for propagating modes, 
±i, for evanescent modes, 



(B17) 



while Cj is generally complex for complex modes, and 
zero for modes with /3j = 0. This is done such that, if 



possible, hfc can be taken to be real. Then (B16) implies 
power orthogonality for these modes: 



hi ■ zdS 



CjS 



•jk- 



(B18) 
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